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ABSTRACT 

Genetic variation in FFAR1 modulates insulin secretion dependent on non-esterified fatty acid (NEFA) concentrations. We previously demon- 
strated lower insulin secretion in minor allele carriers of PPARG Pro12Ala in high-NEFA environment, but the mode of action could not been 
revealed. We tested if this effect is mediated by FFAR1 in humans. Subjects with increased risk of diabetes who underwent oral glucose 
tolerance tests were genotyped for 7 tagging SNPs in FFAR1 and PPARG Pro12Ala. The FFAR1 SNPs rs1 2462800 and rs1 0422744 demon- 
strated interactions with PPARG on insulin secretion. FFAR1 rs1 2462800 (p = 0.0006) and rs1 0422744 (p = 0.001) were associated with 
reduced insulin secretion in participants concomitantly carrying the PPARG minor allele and having high fasting FFA. These results suggest that 
the minor allele of the PPARG SNP exposes its carriers to modulatory effects of FFAR1 on insulin secretion. This subphenotype may define 
altered responsiveness to FFAR1 -agonists, and should be investigated in further studies. 

© 2014 The Authors. Published by Elsevier GmbH. This is an open access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/3.0/). 
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1. INTRODUCTION 

In addition to glucose, which is the foremost stimulator of insulin 
secretion in beta-cells, other nutrients such as free fatty acids (also 
called non-esterified fatty acids, NEFA) contribute to glucose- 
stimulated insulin secretion (GSIS). Glucose triggers the processes 
leading to insulin secretion by its metabolism within beta-cells [1 ,2]. 
NEFA also act as fuels in beta-cells and they are able to modulate 
insulin secretion by influencing intracellular metabolism [3]. It was not 
until 2003, the discovery of the role of a previously orphaned receptor, 
the G-protein-coupled receptor 40 (GPR40), now called free fatty acid 
receptor 1 (FFAR1), that an additional important coupling mechanism 
between NEFA and insulin secretion was revealed [4]. The acute in- 
crease of insulin secretion after raising NEFA levels was shown to be 
~50% in healthy volunteers [5,6], and the FFAR1 -pathway is thought 
to be involved in at least 50% of NEFA-mediated insulin secretion [7]. 
FFAR1 -agonists have been developed to enhance GSIS, some of which 
could be candidates for a new class of antidiabetic drugs [8]. We 
recently demonstrated that FFAR1 -agonism may be protective against 



beta-cell apoptosis and provided evidence that common variation near 
the FFAR1 gene modulates insulin secretion dependent on free-fatty 
acid levels [9]. 

Recent data from an investigation conducted with beta-cell cultures 
and animals suggest that FFAR1 gene expression is modulated by 
stimulation of peroxisome proliferator-activated receptor gamma 
(PPARG) [10]. Intriguingly, in 2001 we had demonstrated effects of the 
PPARG Pro12Ala single nucleotide polymorphism (SNP) on insulin 
secretion in an elevated NEFA milieu. In this investigation, carriers of 
the minor allele of the Pro1 2Ala polymorphism (rs1 8051 92) had lower 
insulin secretion during hyperglycemic clamp studies conducted with a 
concomitant intravenous lipid infusion, but no difference was seen 
between the genotypes without increasing NEFA [11]. 
Given these data supporting a crosstalk between PPARG and FFAR1 
signaling, we set out to look for evidence of a PPARG x FFAR1 
interaction in humans. For this purpose, we analyzed insulin secretion 
effects of interactions between previously described tagging SNPs in 
FFAR1 [9] and the aforementioned diabetes-related SNP in PPARG, 
also accounting for NEFA levels. 
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2. METHODS 

2.1. Study population 

Subjects analyzed in the current study are participants of the Tue- 
bingen Family Study (TUEF). TUEF is a cross-sectional observational 
study originally designed to characterize the phenotype of increased 
risk for type 2 diabetes. Persons with either a positive family history, or 
prior gestational diabetes or known glucose intolerance or overweight 
from Southern Germany were enrolled. For this analysis, we selected 
1928 participants without incident diabetes who had fasting NEFA and 
full 5-point OGTT measurements for glucose and insulin, and who 
were also genotyped for the PPARG and FFAR1 SNPs. Of the study 
population, 3% had prior gestational diabetes, 43% had a family 
history of diabetes in at least one first degree relative and 40% were 
obese according to the WHO definition (BMI > 30 kg/m 2 ) (Table 1). 

2.2. OGTT and laboratory measurements 

Following an overnight fast, all participants ingested a 75 g glucose 
dose at 8 am (OGTT). Plasma glucose, insulin, C-peptide and NEFA 
concentrations were determined after 0, 30, 60, 90 and 120 min. A 
bedside glucose analyzer (glucose-oxidase method, Yellow Springs 
Instruments, Yellow Springs, OH, USA) was used to determine plasma 
glucose. Plasma insulin and NEFA were measured with commercial 
chemiluminescence assays for ADVIA Centaur (Siemens Healthcare 
Diagnostics, Eschborn, Germany). 

2.3. Selection of tagging SNPs and genotyping 

Selection of the tagging SNPs in FFAR1 was performed as described 
earlier [9]. Genotyping was carried out on the MassARFlAY platform 
from Sequenom (Sequenom, San Diego, CA, USA). Initially 9 SNPs 
were selected as tagging SNPs in FFAR1 covering all other common 
SNPs within the locus with an i 2 > 0.8 (1 00% coverage), but genotype 
calling was unsuccessful for rs1 24591 38 and Hardy-Weinberg 
equilibrium was rejected for rs1 0423648 (p < 0.05), therefore 7 SNPs 
went into analysis. 

2.4. Calculations and statistics 

Insulin sensitivity was assessed from 5-point OGTT glucose and insulin 
values using the Matsuda index (ISI), calculated as 
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Table 1 — Demographic and metabolic characteristics of the investigated 
population («= 1928). 



Trait 


Median 


IQR 


Age 


39 


(29, 50) 


BMI (kg/m 2 ) 


28 


(24, 35) 


Glucose 0' (fasting, mmol/l) 


5.1 


(4.8, 5.4) 


Glucose 120' (post-challenge, mmol/l) 


6.2 


(5.2, 7.3) 


HbA1c(%) 


5.4 


(5.1, 5.7) 


Insulin 0' (fasting, pmol/l) 


51.0 


(34.0, 87.0) 


Insulin sensitivity index (AU) 


12.4 


(7.3, 20.5) 


Insulinogenic index (AU) 


124 


(76, 200) 


Non-esterified fatty acids (fasting, u.mol/1) 


561 


(424, 723) 


Triglycerides (fasting, mg/dl) 


100 


(71, 147) 


Cholesterol (fasting, mg/dl) 


190 


(167, 215) 


LDL (fasting, mg/dl) 


116 


(96, 139) 


HDL (fasting, mg/dl) 


52 


(44, 62) 
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where g and /' denote glucose and insulin levels, respectively, at the 
given time-point (minutes) of OGTT [12]. Insulin secretion was 
assessed by the insulinogenic index (IGI) which was originally pro- 
posed by Seltzer et al. in 1967 [13], and has been since then validated 
several times as a reliable measure of insulin secretion utilizing insulin 
and glucose measurements during the first 30 min of an OGTT [1 4,15]. 
IGI can be calculated as (/ 30 - /b)/(£feo - go)- The disposition index (Dl), 
another measure of insulin secretion which also accounts for insulin 
sensitivity, was calculated as IGI x ISI. 

Variables with skewed distributions were log e -transformed prior to linear 
regression analyses. To increase statistical power in the interaction an- 
alyses, genotypes were coded according to the dominant inheritance 
model. Given the 7 investigated SNPs, to reduce the risk of false positive 
findings with multiple testing, the limit of statistical significance was 
defined at 0.05/7 = 0.007, and a p between 0.05 and 0.007 was termed 
nominally significant. Sex, age and insulin sensitivity were used as 
covariates of IGI. Sex and age were used as covariates of Dl. Interaction of 
the FFAR1 SNPs with the Pfl4/?GSNPand with NEFA was tested adding all 
combinations of possible interaction terms (FFAR1 x NEFA, 
PPARG x NEFA, FFAR1 x PPARG, FFAR1 x PPARG x NEFA) to the 
models. Effect estimates are given as standardized beta (/3). All statistical 
analyses were conducted with JMP1 1 (SAS). 

3. RESULTS 

As expected, PPARG genotypes alone did not impact insulin secretion 
in the study population (a table on the genotype-dependent distribution 
of key cohort phenotypes is provided as Supplementary Table 1 ). To 
test the hypothesis that there is an interaction between genetic vari- 
ation in PPARG and FFAR1 on insulin secretion, we analyzed 
gene x gene interactions between rs1 8051 92 in PPARG and 7 tagging 
SNPs of known frequent variants in FFAR1. Since the modulatory effect 
of PPARG on insulin secretion manifested only in a high-NEFA envi- 
ronment [11] and FFAR1 effects on insulin secretion were seen only in 
interaction with fasting NEFA [9], we additionally adjusted for fasting 
NEFA and added interaction terms of fasting NEFA with the FFAR1 
SNPs and the PPARG SNP to the models. 
We adjusted insulin secretion as determined by the insulinogenic index 
for sex, age, insulin sensitivity, fasting NEFA, PPARG Pro1 2Ala geno- 
type, FFAR1 genotypes and 4 interaction terms: between FFAR1 and 
PPARG genotypes, between PPARG genotype and NEFA level, between 
FFAR1 genotype and NEFA and, finally, a 3-way interaction term with 
PPARG x NEFA x FFAR1-SNP. We found a significant gene x gene 
interaction between the FFAR1 SNP rs1 0422744 and the PPARG 
Pro12Ala SNP {p = 0.005) with a concomitant nominally significant 
interaction between NEFA level and the PPARG SNP (p = 0.03). 
Additionally, a nominally significant interaction was found between 
rs1 2462800 and PPARG {p = 0.01) with a concomitant 
PPARG x NEFA interaction (p = 0.03). The 3-way interaction term did 
not reach statistical significance, although the p-value of 0.0592 
suggested a trend for nominal significance in the case of rs1 2462800. 
Detailed results of the interaction tests are provided in Supplementary 
Table 2. An additional adjustment for obesity by adding BMI to the 
interaction models did not relevantly change the results (data not 
shown). Alternatively using Dl as outcome variable, the gene x gene 
interaction term had a p-value of 0.01 and 0.001 for rs1 2462800 and 
rs1 0422744, respectively. 

Next, we conducted subgroup analyses with stratification along 
interaction variables, namely the PPARG SNP and fasting NEFA. The 
study population was thus stratified into high and low NEFA (along the 
median NEFA, 561 Limol/L), and PPARG minor allele (Ala) carriers. This 
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Figure 1: Insulin secretion per rs1 2462800 (A) and rs10422744 (B) genotype in FFAR1 in interaction with the PPARG Pro1 2Ala genotype and fasting free fatty acid levels. The bar 
charts show insulin secretion as residuals of the insulinogenic index adjusted for sex, age and insulin sensitivity (arbitrary units). Dark shaded bar charts represent homozygous 
carriers of the major FFAR1 alleles (G/G for rs1 2462800 and C/C for rs1 0422744), while light shaded bar charts represent hetero- and homozygous carriers of the minor alleles (X/ 
A for rs1 2462800 and X/T for rs1 0422744). The sets of bar charts in the left sided quadrants (I, III) show insulin secretion in the subgroup with the homozygous major PPARG allele 
(Pro/Pro), while the sets of bar charts in the right sided quadrants (II, IV) show insulin secretion in the subgroup that hetero- or homozygously carries the minor PPARG allele (X/Ala). 
The upper (I, II) and lower (III, IV) bar chart sets represent the subgroups with low and high fasting non-esterified fatty acid (NEFA) levels, respectively. This stratification was 
performed along the median NEFA value (561 pmol/l) of the total study population. 
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stratification resulted in 4 subgroups (I: low NEFA and PPARG Pro 
homozygotes, n = 71 8, II: low NEFA and PPARG Ala carriers, n = 237, 
III: high NEFA and PPARG Pro homozygotes, n = 744, IV: high NEFA 
and PPARG Ala carriers, n = 229). An effect of the FFAR1 genotype on 
insulin secretion was only evident in the subgroup IV comprising 
participants with at least one PPARG Ala allele and high fasting NEFA 
(Figure 1A and B bottom right diagrams [diagrams IV]). In this sub- 
population, carriers of the minor allele of rs1 0422744 and rs1 2462800 
in FFAR1 had lower insulin secretion (/? = -0.18, p = 0.001 and 
/5 = -0.19, p= 0.0006, respectively) than homozygotes of the major 
allele, after controlling for sex, age and insulin sensitivity. 

4. DISCUSSION 

In this work, we provide evidence for a biologically relevant crosstalk 
between PPARG and FFAR1 signaling in the regulation of insulin 
secretion in humans. Importantly, this crosstalk was only manifest 
when fasting NEFA levels were high. 

The Pro12Ala variant in PPARG has been one of the first candidate 
SNPs for type 2 diabetes [16,17]. The rare allele of this variant is 
associated with a 25% reduced risk for the disease [1 8]. Given the role 
of PPARG as both a receptor for NEFA and itself a regulator of fat 
metabolism, our group had earlier investigated the hypothesis that 
PPARG-mediated alterations in fatty acid signaling could have an 
impact on insulin secretion. In those hyperglycemic clamp studies, we 
had compared insulin secretion between Pro/Pro and X/Ala carriers of 
the PPARG variant in a control and a high-NEFA setting. Only the clamp 
condition with a concomitant infusion of Intralipid solution, performed 
to raise NEFA levels, had unmasked striking differences in insulin 
secretion between PPARG genotypes [11]. However, no explanation 
could be given for the mechanism of action by that time. Meanwhile, 
the NEFA receptor FFAR1 has been established as an important 
stimulator of fatty-acid mediated insulin secretion [4] and genetic 
variation in FFAR1 has been shown to influence beta-cell function 
[9,19]. Therefore, FFAR1 seemed to be a plausible link between PPARG 
and insulin secretion. By analyzing the previously described 7 FFAR1 
tagging SNPs, we now found 2 SNPs, which exert a NEFA and PPARG- 
dependent effect on insulin secretion. 

The SNPs rs1 2462800 and rs1 0422744 are located 0.8 kb apart in an 
intergenic regulatory area between the FFAR1 (GPR40) and FFAR3 
(GPR41) genes, 3.5 and 3.8 kb from the 3' end of the single FFAR1 
exon. These SNPs are more distal from the gene than the previously 
described FFAR1 SNP rs1 57361 1 which directly interacts with fasting 
NEFA in association with insulin secretion [9]. 
The molecular mechanisms underlying the PPARG-FFAR1 interaction 
are still speculative. Although the literature is controversial, there is a 
wealth of data from basic science and clinical studies involving the use 
of PPARG agonists (thiazolidinediones) indicating that PPARG activation 
may have a positive impact on beta-cell function and beta-cell survival 
(reviewed by [20]). The PPARG-FFAR1 interaction, especially in light 
of the effect modification by NEFA levels, could possibly explain this 
still poorly understood link between PPARG action and beta-cell 
function. A recent study demonstrated that PPARG overexpression 
leads to an increased expression of FFAR1 in rat INS-1 cells and 
primary rat islets [10]. A possible interaction scenario could thus be a 
defective PPARG-induced transcriptional or translational activity of the 
FFAR1 gene in the presence of the FFAR1 minor allele variants, which 
would consecutively result in a limited compensatory potential to in- 
crease insulin secretion when fasting NEFA are elevated. The funda- 
mental human phenotypic alterations associated with the Pro12Ala 
variant seem to involve lipolysis and its insulin-dependent inhibition 
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[21-23]. Since FFAR1 is activated by medium and long-chained fatty 
acids [24], perturbations in NEFA levels or their composition as a 
consequence of altered PPARG action could modulate insulin secretion 
through FFAR1 . 

In summary, PPARG and FFAR1 are connected in several ways, and the 
presented interaction is mechanistically reasonable. Additional in- 
vestigations are required to elucidate the exact details of the underlying 
physiology. However, our findings are of immanent importance for the 
treatment of type 2 diabetes, because the described genotypic inter- 
action between PPARG and FFAR1 could have bearing on the phar- 
macologic efficacy of FFAR1 agonists. Although a phase 3 clinical trial 
of the most promising FFAR1 -agonist has been recently halted, several 
of such agents are still under development [8,25]. We infer that the 
variants identified in the current study, under the constellation of the 
minor PPARG Pro12Ala allele, could have relevant impact on the 
effectiveness of FFAR1 agonists. Identification of variants which are 
candidates for pharmacogenetic interactions could pave the way for 
investigating genotype-based therapeutic responses within clinical 
studies and on the long run in clinical practice. This would eventually 
lead to an improved selection of individuals who would benefit from 
specific therapies [26]. 
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